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HS 2237+8154: On the onset of mass transfer 
or entering the period gap?* 
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Abstract. Wc report follow-up observations of a new white dwarf/red dwarf binary HS 2237+8154, identified as 
a blue variable star from the Hamburg Quasar Survey. Ellipsoidal modulation observed in the R band as well 
as the Ha radial velocity variations measured from time-resolved spectroscopy determine the orbital period to 
be Porb = 178.10 ± 0.08 min. The optical spectrum of HS 2237-1-8154 is well described by a combination of a 
Tcff = 11500 ± 1500 K white dwarf (assuming logg = 8) and a dM3.5±0.5 secondary star. The distance implied 
from the flux scaling factors of both stellar components is d = 105 ± 25 pc. Combining the constraints obtained 
from the radial velocity of the secondary and from the ellipsoidal modulation, we derive a binary inclination of 
z ~ 50° - 70° and stellar masses of Mwd = 0.47 - 0.67 and Msec = 0.2 - 0.4 M©. AU observations imply 
that the secondary star must be nearly Roche-lobe filling. Consequently, HS 2237-1-8154 may be either a pre- 
cataclysmic variable close to the start of mass transfer, or - considering its orbital period - a cataclysmic variable 
that terminated mass transfer and entered the period gap, or a hibernating nova. 
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1. Introduction 

The currently known population of cataclysmic variables 
(CVs) comprises more than 1000 system s, with period 
measurements available for :^ 500 of them jPownes et al.l 
I2n0lt iKube et al.ll20n,'4 bitter fc Kolbl l2nn.'^ . Compared 
to these numbers, on ly very few CV pro genitors have 
been identified so far. ISchreiber fc Gansicka (,2003,) have 
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recently analysed 30 well-observed pre- CVs and found 
that the sample of known pre-CVs is strongly biased to- 
wards systems containing hot (young) white dwarfs and 
late type (low-mass) companion stars. Testing and im- 
proving our understanding of the evolution of pre-CVs 
through the common envelope phase and of the angular 
momentum loss mechisms that subsequently brings them 
into a semidetached configuration will enormously benefit 
from a larger and unbiased sample of such stars. 



We have identified a detached white dwarf/red dwarf 
binary (Fig.nj in the course of our studies of the 
stellar content of the 



Gansicke ct al. 


2000; 


Nogami ct al. 2000; 


Szkodv et alJ 


2001 : ICansicke 


et al.ll2002l: lArauio-Betancor et al.ll2003l). 



Our follow-up observations show that HS 2237-1-8154 may 
belong to the large population of " c old fc old" pre-CVs 
predicted by ISchreiber fc Gansickd l)2003(l . however, it 
may as well be a CV in a peculiar phase of its evolution. 
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2. Observations and Data Reduction 

2.1. Spectroscopy 

A single identification spectrum of HS 2237+8154 was ob- 
tained in June 1991 using the 3.5m telescope at Calar 
Alto equipped with the TWIN spectrograph (Tabled]). 
Gratings with 144 and 160 A/mm were used in the blue 
and red channel respectively providing a resolution of 
about 6 A. Standard data reduction with flat fielding, bias 
subtraction as well as cosmic-ray filtering was performed 
in Kiel using IDAS, a long-slit spectroscopy package writ- 
ten by G. Jonas (University Kiel). The TWIN spectrum 
(Fig.lSJ of HS 2237+8154 clearly reveals the composite na- 
ture of the object, with prominent TiO absorption bands 
in the red part of the spectrum, hallmarks of a late 
type star, and broad Balmer absorption lines in the blue 
part of the spectrum, characteristic of the high-gravity 
atmosphere of a white dwarf. Narrow emission of Ha 
and H/3 is detected, which is most likely associated with 
chromospher ic activity on the late-type companion (e.g. 
iBleach et al l2000'). 

Additional time-resolved spectroscopy of 
HS 2237+8154 was obtained in August 2002 with 
the Intermediate Dispersion Spectrograph (IDS) on the 
Isaac Newton Telescope, totalling 30 spectra spread out 
over one week (TableQ]). The IDS was equipped with 
the R632V grating and the 2048 x 4100 pixel EEVlOa 
detector. Using a slit width of 1.5"the setup provided an 
unvignetted wavelength range of ~ 4400 A — 6800 A and a 
spectral resolution of ~ 2.3 A. Copper-argon wavelength 
calibrations (arcs) were obtained at the beginning and 
end of each observation block of HS 2237+8154. The 
IDS data were bias-subtracted and flat-fielded in a 
standard manner using the Figaro package wit hin the 
Starlink software collection. Optimum extraction l(Hornel 
Il986l') and sky line subtraction was carried out using 
Tom Marsh's (^989) Pamela package. The dispersion 
relation was established fitting the arc line positions 
with a fifth-order polynomial. The rms of the wavelength 
solution was smaller than 0.1 A for all spectra. Possible 
drift of the wavelength-to-pixel scale was accounted for 
by interpolation between the arcs. 

2.2. Photometry 

We have obtained differential i?-band photometry of 
HS 2237+8154 during five nights in September/October 
2002 at the 1.2 m Kryoneri telescope using a SI-502 
516 X 516 CCD camera (Table^. Bias and dark-current 
subtraction as well as flat-fielding of all CCD images was 
carried out in a standard fashion within MIDAS. Because of 
the relatively poor tracking none of the standard photom- 
etry packages could be used to produce light curves from 
the CCD data without excessive human intervention, and 
we therefore set up the following reduction pipeline. (1) 
We extracted detector coordinates and instrumental mag- 
nitudes for all objects in each image using the sextractor 
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Fig. 1. Finding chart (10' x 10') for HS 2237+8154 ob- 
tained from the Digitized Sky Survey. The coordinates 
of the star are a(2000) = 22''37™15.8^ and (5(2000) = 
+82°10'27.5". The comparison stars used in the analysis 
of the Kryoneri data are marked 'CI' (USNO-A2.0 1650- 
02587295 ) and 'C2' (USNO-A2.0 1650-02586211). 

Table 1. Log of the observations. 



Date 


UT Time 


Data 


Exp.(s) 


Num. Obs 


Spectroscopy 




ii~6A 






1991 Jun 24 


03:30-03:45 


900 


1 


2002 Aug 28 


00:32-01:55 


i? ~ 1.6 A 


600 


9 


2002 Aug 30 


01:35-02:58 


i? ~ 1.6 A 


600 


8 


2002 Sep 01 


01:51-02:43 


i? ~ 1.6 A 


600 


6 


2002 Sep 03 


23:35-01:31 


i? ~ 1.6 A 


600 


7 


Photometry 










2002 Sep 19 


18:31 - 23:40 


R 


10 


961 


2002 Sep 20 


17:42 - 22:03 


R 


10 


960 


2002 Oct 08 


17:16 - 22:17 


R 


10 


982 


2002 Oct 17 


17:38 - 21:52 


R 


30 


415 


2002 Oct 18 


17:00 - 23:18 


R 


45 


417 



I Bertin fc Arnoutl Il996l) in the "classic" aperture pho- 
tometry mode. In order to account for variable seeing, this 
step is done as a two-stage process: sextractor is run in 
a first pass using a relatively large aperture diameter (10 
pixel), the median seeing is computed from the resulting 
object list, and then sextractor is run a second time, us- 
ing an aperture diameter of 1.5 times the median seeing, 
which proved to give the best signal-to-noise ratio. (2) One 
CCD image is designated to be the "reference frame" , and 
the object lists derived from all other images are matched 
against that of the reference frame. Considering that most 
of our CCD images of HS 2237+8154 have only three ob- 
jects in common, and that no significant distortion/scaling 
occurs between two different images, we decided to imple- 
ment the purpose-build matching algorithm smatch. This 
program computes the distances between all possible pairs 
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Fig. 2. The Calar Alto TWIN idenfication spectrum 
of HS 2237+8154. Plotted in gray is a two-component 
fit consisting of a white dwarf model {Tctt — 12 000 K, 
logg = 8.0) and a dM3.5 template (G1273). 



of stars in a given image, and compares them with the dis- 
tances of all possible pairs of stars found in the reference 
image. A "match" is found if the distances between two 
stars agree within less than a given uncertainty e, with 
e < 3 pixel in the case of the HS 2237+8154 data, smatch 
proved to be robust for small numbers of stars (5 — 6) and 
provided a better match-rate than more complex codes 
such as IRAF's xyxymatch. (3) As a consequence of the 
small CCD detector in use at the Kryoneri telescope (with 
a field of view of 2.5' x 2.5') and the poor tracking, special 
attention has to be given to the choice of the comparison 
stars, as some of them may move off the detector during 
the observations. We therefore assign up to four different 
comparison stars, and compute the differential magnitudes 
of the target star relative to the brightest comparison vis- 
ible on a given image. 

i?-band magnitudes of HS 2237+8154 were derived rel- 
ative to USNO-A2.0 1650-02587295 (R = 14.5) for the 
observations obtained in September 2002, and relative to 
USNO-A2.0 1650-02586211 {R = 15.6) for the observa- 
tions obtained in October 2002 (see Fig^. During all five 
nights, HS 2237+8154 was found at a mean magnitude of 
R = 15.73±0.02. However, it is important to keep in mind 
that the systematic error on this R magnitude is of the 
order ^ 0.2 mag due to the limited photometric accuracy 
of the USNO-A2.0 catalogue (which hsts HS 2237+8154 
with B = 16.1 and R = 15.9). For completeness we note 
that B = 16.1 has been derived from the HQS observa- 
tions of HS 2237+8154. 

The light curves of HS 2237+8154 (Fig.Ej) display a 
quasi-sinusoidal variation of the i?-band magnitude with 
a peak-to-peak amplitude of ^ 0.1 mag and a period of 
~ 90min during all five nights. 




2637.3 2537.35 2537.4 2537.45 2537.5 




253B.25 2538.3 2538.35 2538.4 2538.45 2538.5 




2558.25 2556.3 2556.35 2556.4 2556.45 



2002-10-17_ 



2565.25 2565.3 2565.35 2565.4 2565.45 



2002- la- 18 



2566.25 2566.3 2568.35 2586.4 2568.45 

JD - 24 500 000 

Fig. 3. Differential CCD i?-band photometry obtained at 
the Kryoneri 1.2 m telescope. The R magnitudes have been 
computed relative to the comparison stars 'CI' (09/19 and 
09/20) and 'C2' (10/08 to 10/18). The data from the first 
three nights have a lower signal-to-noise ratio because of 
the relatively short exposure times (10s). 



3. The orbital period of HS 2237+8154 

3.1. Radial velocities 

We have measured the radial velocity variation of the 
narrow Ha emission line in the INT spectra using 
Gaussian fits. Care was taken to avoid a contamina- 
tion of the measurement by the broad M-dwarf con- 
tinuum peak located close to Ha. The radial velocities 
were t hen sub.jected to the analysis- of-variance algorithm 
fAOV ISchwarzenberg-Czernvlll98^ . as implemented in 
the MIDAS context TSA. The resulting periodogram con- 
tains a number of peaks in the range ~ 2 — 10 d~^, with 
the strongest signal at 8.085 d~^ (Fig.EJ. We have cre- 
ated a faked set of radial velocities, computing a sine 
wave with a frequency of 8.085 d~^and adopting the tem- 
poral sampling as defined by the observed spectra. The 
errors of the observed radial velocity measurements were 
used to randomly offset the faked radial velocities from 
the value of the sine wave. The AOV analysis of this 
faked data set is also shown in Fig.^ Whereas the spec- 
troscopy strongly suggests an orbital frequency (period) 
of 8.085 d^^ (178min), it does not permit a secure identi- 
fication. 

3.2. Photometry 

The AOV periodogram computed from our photometric 
data of HS 2237+8154 is shown in Fig.|Sl The strongest 
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Fig. 4. Main panel: analysis-of-variance periodogram of 
the radial velocities measured from the Ha emission line. 
Shown on top is the AOV periodogram of a faked set of 
radial velocities (see text for details). 
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Fig. 5. Main panel: analysis-of-variance periodogram 
computed from the five nights of differential photometry 
shown in Fig.El On top: the AOV periodogram computed 
from a sine wave with / = 2 x /orb reproduces all struc- 
tures detected in the data. 



peak is found at / ~ 16 d^^ [P ~ 90min), which cor- 
responds to the period of the quasi-sinusoidal modula- 
tion apparent in the light curves (Fig.UJ. A second struc- 
ture in the periodogram is centred at half that frequency, 
~ 8d~^ {P ~ 180 min). Considering that the i?-band flux 
of HS2237-f8154 is almost entirely dominated by emis- 
sion from the secondary star, the observed quasi-sinusoidal 
variability is readily explained by either ellipsoidal mod- 
ulation, in which case Porb — 180 min, or by reflection 
effect of the secondary irradiated by the white dwarf, in 
which case Porb — 90 min. The absence of a signal in the 
~ 16 d~^ range of the periodogram obtained from the ra- 
dial Ha velocities (Fig.0 clearly rules out the second op- 
tion. We therefore identify the ~ 8d~^ signal, which is 
detected consistently both in the photometric variability 
as well as in the radial velocity variation, as the orbital 
period of HS 2237-f 8154. 

Combining the photometric and spectroscopic obser- 
vations, we define the following ephemeris 



HJD 2452514.4112(1) + 0.12368(6) x E 



(1) 



were phase zero is the time of the blue-to-red crossing 
of the Ha radial velocities. The most accurate measure- 
ment of the orbital period, Porb = 178.10 ± 0.08 min, is 
derived from the photometry which covers a longer time 
span compared to the spectroscopy. 

FiguresEl and \7\ show the Kryoneri photometry and 
the INT radial velocities folded with this ephemeris. The 
relative phasing of the photometry and spectroscopy is 
consistent with our interpretation of the photometric vari- 
ability as ellipsoidal modulation: photometric minima oc- 
cur during the zero-crossing of the radial velocity curve 
(i.e. when the front/back side of the secondary star is 
seen), photometric maxima occur during the radial ve- 
locity maxima (i.e. when the secondary star is seen from 



the side). A trailed spectrogram of the phase-folded IDS 
spectra around Ha is shown in Fig.|51 

4. The stellar components in HS 2237+8154 

4.1. Spectral fit 

We have fitted the identification spectrum of 
HS 2237+8154 (Fig.EJ with a two-component model, 
consisting of synthetic whit e dwarf spectr a calc ulated 
with the code described by ICansicke et al.l lll995l). and 
observed M-dwarf templates from iBeuermann et alJ 
(1998). We fixed the surface gravity of the white dwarf 
models to log 5 = 8, corresponding to a canonical white 
dwarf mass of ^ 0.6 Mq, and allowed for temperatures 
in the range 8000 - 25 000 K. The spectral hbrary of 
M-dwarfs extended from MO. 5 to M9. The closest match 
of the observed spectrum was achieved for a white dwarf 
with Tcff = 11 500 ± 1500 K and a dM3.5±0.5 companion 
(Fig. 121). The determination of the secondarys spectral 
type is relatively robust, as the late- type star completely 
dominates the emission of the system in the red part of 
the spectrum. The flux from the white dwarf exceeds that 
of the companion only for A ^ 5000 A, and consequently 
its temperature is subject to a relatively large uncertainty. 
The scaling factors between the model/template spectra 
and the observed spectrum of HS 2237+8154 allow an 
estimate of its distanc e. Fro m the white dwarf (assuming 
a Hamada-Salpeter lll96ll) mass-radius relation), we 
obtained d — 115 ± 15 pc. Obviously, the distance would 
be lower (higher) if the white dwarf were more (less) 
massive than 0.6 Mq. Assuming that the companion 
star fills its Roche volume (see the discussion in Sect.O, 
the scaling factor of the M-dwarf template implies 
(i = 95 ± 15 pc. The agreement between the two values 
is satisfying, and considering the overall uncertainties 



B.T. Gansicke et al.; HS 2237+8154: On the onset of mass transfer or entering the period gap? 



5 



involved in the spectral fit d = 105 ± 25 pc appears 
to be a conservative estimate, well within the reach of 
groun d-based astrometry programs (e.g. iThorstenseiJ 
120031) . 



4.2. The radial velocity of the secondary 

Fitting a sine curve to the phase-folded radial velocities 
(Fig.[7I) gives a semi-amplitude Kscc — 236 ±8 km s^^ and 
a systemic velocity 7 = —40 ± 6kms~^. Adopting this 
value of Kscc implies a mass function 
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Possible (Mwd, Msec) Combinations that satisfy this mass 
function are plotted for a variety of inclinations in Fig.|^ 
Inclinations i ^ 75° are excluded from the absence of 
eclipses (assuming that the secondary is nearly filling 
its Roche lobe, see Sect.|5)l. The typical mass range for 
field M-dwar fs with a spectral type dM3.5ib0.5 is Ms^.r. — 
0.2 -- 0.4M,7>llKirkDatrick et al.lll99ll: iBaraffe fc Chabrieil 
llflflfitlDelfosse et a,1.lll99(£ 

4.3. Ellipsoidal modulation 

The observed i?-band light curves show an ellip- 
soidal modulation with a peak-to-peak amplitude of ~ 
0.085 mag. No evidence for a reflection effect is detected - 
not too big a suprise considering that the white dwarf is 
rather cold. Intriguing is the fact that the observed max- 
ima in the i?-band light curve are not equally bright. A 
similar effect has been obse rved in the i?-band light curve 
of the pre-CV BPM 71214 feawk a fc Yenned l2003|) . who 
suggested an inhomogenous distribution of starspots as 
possible explanation. In fact, rapid rotation due to tidal 
locking may result in enhanced activity on the secondary 
stars of pre-CVs. 

We have computed ellipsoidal modulation model light 
curves for a grid of (Mwd. Mscc^ i) combinations using the 
analytical formulae of Binnendiik (1974), adopting gravity 
and limb darkening coefficients of y = 1.38 and x = 0.80 
(the detailed choice of these two coefficients does not sig- 
nificantly alter the results described below). More specif- 
ically, we have searched (M^d, -^^soc) combinations that 
reproduce the observed amplitude of the ellipsoidal mod- 
ulation for a fixed value of i. In first order approximation 
this implies a constant mass ratio q — Afgcc/Afwd- The 
resulting constraints on (M^d, Msec) are shown in Fig.|51 
for i = 50° and i = 70°. 

Combining the constraints from observed amplitudes 
of the ellipsoidal modulation and the Ha radial velocity 
variation provides a lower limit on the Roche-lobe fill- 
ing factor of the secondary. Assuming the highest incli- 
nation and highest mass ratio that are consistent with the 
mass function f{Mscc) {i ^ 70°, q = M^cc/Mwd ^ 0.7) 
implies i?sGc ^ 0.75i?i,i (with Rl^ the equivalent radius 
of a Roche-lobe filling secondary). Decreasing i and/or q 



Fig. 8. Trailed spectrogram of the Ha emission line in the 
IDS spectra of HS 2237-1-8154. The 30 individual spectra 
have been normalised and averaged into 20 phase bins. 



:) requires i?soc to be closer 
i 60°, and i?scc — Rli for 



within the constraints of /(Msec; 
to Rl^ (Rscc ^ 0.9Rli for i : 
i ~ 50°). 

If the secondary star in HS 2237-1-8154 has a mass that 
is typical for its spectral type, and nearly fills its Roche 
volume, then the combination of all constraints shown in 
Fig. 1^1 results in conservative estimates of the binary in- 



clination and of the component masses: 



50° - 70° 



Ms, 



0.2 - 0.4 Mq, and M, 



0.47- 0.67 M, 



0- 



5. The evolutionary state of HS 2237+8154 

Based on its composite spectral appearence and orbital 
period HS 2237-f 8154 qualifies as a bona-fide pre-CV. 
However, the spectral type of the secondary is con- 
spicously cl ose to that found in C V with orbital periods 
around 3 h l)Beuermann et al.ll998(l . Could HS 2237-1-8154 
be a genuine CV caught in a state of low accretion? We be- 
lieve that this is rather unlikely for the following reasons. 
(1) The system has been found near S ~ i? ~ 16 dur- 
ing six epochs of observations. (2) A number of quiescent 
dwarf novae are known to reveal the photospheric emission 
of their stellar components during quiescence. However 
their optical spectra are always significantly contaminated 
by the accretion disc, in the form of strong Balmer lines 
and some m ore or less noticeabl e contribution to the con- 
tinuum fe.g. lSzkodv et al"ll2000(l . No emission from a qui- 
escent disc is observed in HS 2237-1-8154. (3) Many of 
the strongly magnetic cataclysmic variables (polars) ex- 
hibit deep low states during which their spectra are com- 
pletely dominated by the stellar components (e.g. AM Her 
.Gansicke et al . 1995; EFEri Bcuermann et al. 2000) . The 
non-detection of Zeeman-splitting in the H/J to Hi5 absorp- 
tion lines limits a hypothetical field of the white dwarf in 
HS 2237+8154 to B <, 5MG , well below the lowest field 
measured in polars. 

Whereas the long-term behaviour and the spectral 
characteristics argue against a genuine CV nature of 
HS 2237+8154 the spectral type of its secondary combined 
with the observed constraints from the observed ellipsoidal 
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Fig. 6. The Kryoneri photometry average-subtracted and Fig. 7. Radial velocties measured from the Ha emission 

folded using the ephemeris Eq.n(gray points) and binned line, folded over the ephemeris given in Eq.^ Plotted as 

into 25 phase bins (black points). Phase zero is defined dashed line is a sine fit to the radial velocity data which 

as the blue-to-red crossing of the Ha radial velocity (see gives a semi-amplitude Kscc = 236 ± 8kms~^ and a sys- 



Fig.EI). 



temic velocity 7 ~ —40 ± 6 km s 



modulation and radial velocity variation (Fig.^ clearly 
show that the secondary star must be nearly filling its 
Roche volume. 

With a temperature of Tcff = 11 500 K the white dwarf 
in HS 2237-1-8154 is unusually cold compared to the other 
known pre-CVs and with an implied cooling age of ^ 400 
million years (fo r Mwh — 0.6 Mrr^) it is one of the old- 
est pre-CVs ('see ISchreiber fc Gansick3l2003l) . Taking the 
available results at face value, we consider three possible 
evolutionary scenarios for HS 2237-1-8154. (1) It is a pre- 
CV close to the onset of mass transfer. Ironically enough, 
it will be "born" close to the upper edge of the period 
gap. (2) It is a well-behaved CV - in the context of the 
disrupted magnetic braking model of CV evolution - that 
turned off mass transfer and entered the period gap. (3) It 
is a CV that is currently hibernating in an extended low 
state, e. g. as a consequence of a moderately recent nova 
eruption IShara et 

It is interesting to note that two other pre-CVs 
have recently been suggested to be hibernating novae 
as they contai n nearly Roche-lobe fi lling secondaries: 
EC 13471-1258 J O'Donoghuc et al.' 20031 and BPM 71214 
( Kawka fc Vennea.2003,,) . For these two systems, our op- 
tion (2) is not available, as they have Porb = 217 min and 
Porb — 290 min, respectively. 

6. Conclusions 

We have identified a detached white/red dwarf binary with 
an orbital period of Porb ~ 178.10±0.08min. Our spectro- 
scopic and photometric observations constrain the binary 
inclination to be z c± 50° — 70° and the stellar masses 
to be Mwd = 0.47 - 0.67 Mq and M,oc = 0.2 - 0.4 Mq. 
Modelhng the optical spectrum of HS 2237-1-8154 we find 
Toff = 11 500 ± 1500 K for the white dwarf and a spectral 
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Fig. 9. Constraints on the binary inclination i and the 
stellar masses. Msec and Mwd- The gray shaded area is 
the mass range typically found in isolated dM3.5±0.5 
stars. The solid lines indicate (M^d, -Msec, «) combinations 
consistent with the measured mass function/ (Msec) = 
0.17 Mq. Plotted as dashed lines are constraints derived 
from modelling the observed ellipsoidal modulation. 



type dM3.5±0.5 for the secondary. The distance of the 
system implied by the flux scaling factors of both stellar 
components is d = 105 ± 25 pc. 

The observations strongly suggest that the secondary 
star in HS 2237-f 8154 is nearly filling its Roche lobe, and 
we conclude that the system is either a pre-CV just about 
to evolve into a semi-detached configuration, a CV that 
evolved into a detached configuration and entered the pe- 
riod gap, or a hibernating nova. All three possible scenar- 
ios are extremely interesting in the context of CV evolu- 
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tion, and we encourage future work to refine the binary 

parameter of HS 2237+8154. 
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